Hydrogen Embrittlement of Aluminum: the Crucial Role of Vacancies 
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We report first-principles calculations which demonstrate that vacancies can combine with hydro- 
gen impurities in bulk aluminum and play a crucial role in the embrittlement of this prototypical 
ductile solid. Our studies of hydrogen- induced vacancy superabundant formation and vacancy clus- 
terization in aluminum lead to the conclusion that a large number of H atoms (up to twelve) can 
be trapped at a single vacancy, which over-compensates the energy cost to form the defect. In the 
presence of trapped H atoms, three nearest-neighbor single vacancies which normally would repel 
each other, aggregate to form a trivacancy on the slip plane of Al, acting as embryos for microvoids 
and cracks and resulting in ductile rupture along the these planes. 
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Hydrogen degradation of the structural properties of 
solids, referred to as embrittlement, is a fundamental 
problem in materials physics. Despite intense studies, 
the definitive mechanism of H embrittlement in metals 
remains poorly understood. Four general mechanisms 
have been proposed: (i) formation of a hydride phase; 
(ii) enhanced local plasticity; (iii) grain boundary weak- 
ening and (iv) blister and bubble formation The un- 
derlying atomic processes and relative importance of the 
four mechanisms remain uncertain, and it is likely that 
a combination of these processes may contribute to em- 
brittlement simultaneously. For these mechanisms to be 
operational, however, a critical local concentration of H 
is required, either to form a hydride phase or to initi- 
ate cracking at microvoids and grain boundaries. One of 
the outstanding problems in the current theories of hy- 
drogen embrittlement is the lack of a comprehensive and 
coherent atomistic mechanism to account for the critical 
H concentrations at crack tips. Moreover, it is widely 
observed that H-cnhanced dislocation mobility is a pre- 
lude to the embrittlement and that the fracture planes 
coincide with the slip plane of the material, which is not 
the typical situation |l|; now ai l these phenomena come 
about still remains a mystery. It is generally believed 
that dislocations are central to H embrittlement phenom- 
ena, and a large body of work has been dedicated to 
elucidate hydrogen-dislocation interaction and its conse- 
quences on embrittlement Vacancies, being ubiq- 
uitously present in solids and having the ability to act as 
impurity traps, could play a central role in the embrittle- 
ment process, but detailed arguments about this role or 
estimates of its relative importance are totally lacking. 

Recent experiments on H-metal systems offer clues on 
the role that vacancies may play in H embrittlement. One 
set of experiments has established that H could induce 
superabundant vacancy formation in a number of met- 
als, such as Pd, Ni, Cr, etc. 0,0. The estimated vacancy 
concentration, Cy, in these systems can reach a value 
as high as 23 at.% ; 3f. A conclusion drawn from these 
experiments is that H atoms, originally at interstitial po- 
sitions in the bulk, are trapped at vacancies in multiple 



numbers with rather high binding energies. It was spec- 
ulated that several (three to six) H atoms can be trapped 
by a single vacancy, with the highest number (six) cor- 
responding to the number of octahedral sites around a 
vacancy in either the fee or the bee lattice Q. Actually, 
we shall show below based on first-principles theoretical 
calculations that in Al, the prototypical simple metal and 
ductile solid, up to twelve H atoms can be trapped at a 
single vacancy site. The consequence of H trapping is 
that the formation energy of a vacancy defect is lowered 
by a significant amount, an energy that we define as the 
H trapping energy. Such reduction in the vacancy for- 
mation energy could result in drastic increase (10 7 fold 
for Fe) of equilibrium vacancy concentrations |5J. The 
superabundant vacancy formation in turn provides more 
trapping sites for H impurities, effectively increasing the 
apparent H solubility in metals by many orders of magni- 
tude. For example, it was observed experimentally that 
about 1000 atomic parts per million (appm) of H atoms 
can enter Al accompanied by vacancy formation at the 
surface under aggressive H charging conditions, which 
should be contrasted with the equilibrium solubility of H 
in Al of about 10 -5 appm at room temperature where 
the experiments were carried out|(|; this is a staggering 
change of eight orders of magnitude in concentration. It 
was futher observed that the H- vacancy defects clustered 
and formed platelets lying on the {111} planes, which di- 
rectly lead to void formation or crack nucleation on the 
{111} cleavage planes 0- 

In order to elucidate the complex nature of H- vacancy 
interaction and to shed light on experimental results, 
we have performed first-principles calculations to exam- 
ine the energetics and electronic structure for the rele- 
vant H- vacancy complexes in Al. Due to the extremely 
low solubility of H in bulk Al, experiments are usu- 
ally difficult and results are dependent on H charging 
conditions; for such systems, first-principles calculations 
are particularly useful to complement experimental ap- 
proaches. Our first-principles calculations are based on 
density functional theory with the VASP implementation 
7] and ultra-soft pseudopotentials ||. The local-density 
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FIG. 1: Schematic representation of the environment of a 
vacncy in Al. (a) The vacancy as a large open circle and 
its 12 nearest neighbors as smaller grey circles, which lie on 
highlighted [100] planes. The cube in dashed lines represents 
the conventional cell of the FCC lattice of side a. A shaded 
tetrahedron with one corner at the vacancy site is also shown, 
and the geometric center of which corresponds to the lowest- 
energy site for a H interstitial atom (shown as black circle) in 
bulk Al. (b) The arrangment of four of the six H2 molecules 
surrounding the vacancy, on a [100] plane, with the first- and 
second-nearest-neighbor Al atoms indicated. The other two 
molecules lie directly above and below the plane of the figure, 
along an axis perpendicular to this plane passing throught 
the vacancy site. In both (a) and (b) the ions are placed at 
the ideal lattice sites, with the atomic relaxations not shown 
explicitly. 

approximation (LDA) is used in all of our calculations, 
with checks based on the generalized gradient approxi- 
mation (GGA) for selected cases. For Al, we find that 
LDA results are consistently closer to experimental val- 
ues than GGA results, so here we will rely mainly on 
LDA numbers to draw physical conclusions. We employ 
a supercell containing 108 atomic sites in a simple cubic 
lattice to model bulk Al, with a 4 x 4 x 4 reciprocal space 
grid in the suprecell Brillouin zone and a plane-wave ki- 
netic energy cutoff of 220 eV for the Al-H system. With 
these parameters, we obtain the formation energy of a 
single vacancy (0.66 eV) and the binding energy for the 
nearest-neighbour (NN) divacancy in pure Al j-0.06 eV) 
in excellent agreement with other theoretical [9j and ex- 
perimental results [13 (Table 1). We note that the NN 
divacancy formation energy is negative, implying that 
it is unstable compared to two isolated single vacancies. 
This counter-intuitive result is due to charge redistribu- 
tion in the neighborhood of the vacancy, which has been 
interpreted as formation of directional covalent/metallic 
bonds that stabilizes the single vacancy configuration 
against the formation of the divacancy 0, El • 

Our main objective is to understand the atomistic 
mechanisms of H- vacancy interaction in Al. First we ad- 
dress the relative site preference of H in bulk Al. To 
this end, we have calculated the total energy of a single 
H atom situated near the vacancy site, or at intersti- 



tial tetrahedral and octahedral bulk sites which are as 
far as possible from the vacancy within the supercell. 
For H atoms, the tetrahedral interstitial site in bulk Al 
is slightly more favorable than the octahedral intersti- 
tial site by 0.07 eV. We find that the H atom prefers to 
occupy the vacancy site over the interstitial tetrahedral 
site in bulk by 0.40 eV. The corresponding experimental 
value is 0.52 eV Q], and theoretical results range from 
0.33 eV to about 1 eV GUI!. The lowest energy posi- 
tion for the H atom in the presence of a vacancy is not 
at the geometric center of the vacancy site, but rather at 
an off-center position close to a tetrahedral site adjacent 
to the vacancy site (see Fig. 1(a)); the energy differ- 
ence between the center and off-center positions is 0.66 
eV. We also find that the H atom is negatively charged, 
consistent with the view that the H impurity can be re- 
garded as a screened H ion in free-electron-like metals 
|l4|. Previous studies based on the jellium model of Al 
have shown that as the jellium conduction electron den- 
sity decreases, the excess charge buildup at the H atom is 
also reduced and the electrons of the H~ ion are less lo- 
calized 0,^3- Therefore, the kinetic energy of the H~ 
electrons is lowered at the vacancy site where the con- 
duction electron density is lower. At the same time, it is 
energetically favorable for the H~ ion to sit off-center of 
the vacancy, to minimize the Coulomb interaction energy 
with the nearby Al ions. 

Having established the stability of a single H atom at 
a single vacancy in Al, the ensuing question is whether 
multiple H atoms, in particular, H2 molecules would be 
stable at this defect. This is an interesting problem on 
its own right, but it is also relevant to H2 bubble forma- 
tion that gives rise to H embrittlement. To examine the 
stability of an H2 molecule at a vacancy site, we com- 
pare the binding energy of the H2 unit at the vacancy 
and in vacuum. The binding energy E b of the H2 unit at 
a vacancy site is calculated as: 

E b = E C (V M + H 2 ) + E C (V M ) - 2E c (y Ai R), (1) 

where E c (Vai + H 2 ) is the cohesive energy of a system 
with an H2 unit at the center of the vacancy, E c (Vai) is 
the cohesive energy of a system with a single vacancy in 
the absence of the H2 unit, and -E c (VaiH) is the cohesive 
energy of a system with a single H atom at the vacancy 
(in the off-center tetrahedral site). Interestingly, we find 
the this binding energy to be +0.06 eV, indicating a weak 
repulsion between the two H atoms in the H2 unit at the 
vacancy site. This is to be compared to the binding en- 
ergy of an H2 molecule in vacuum, which is —6.67 eV. 
The positive binding energy of H2 at the vacancy site 
does not imply that there is no bonding between the two 
H atoms; it simply states that these two H atoms would 
prefer to be trapped at two single vacancy sites individu- 
ally rather than in the same vacancy site as a pair. The 
weakening of the H-H bond at the vacancy site is remark- 
able given the fact that each H atom in the H2 unit in 
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FIG. 2: Trapping energy per H atom in eV as a function of 
the number of H atoms being trapped at a single vacancy site. 
The zero energy corresponds to the energy of a H atom at the 
tetrahedral interstitial site. 



this situation is quite far away (2.6 A) from the nearest 
Al ions. 

We find that the equilibrium interatomic distance be- 
tween the H atoms at the vacancy is 0.83 A, 12% longer 
than the H2 bond length (0.74 A) of the molecule in vac- 
uum; this is due to the partial occupation of antibonding 
states between the H atoms. This can be understood as 
follows: each H atom is associated with a doubly occu- 
pied bound state in the presence of conduction electrons, 
and hence is negatively charged. When the two H atoms 
approach each other, the two bound states split up into 
a bonding and an antibonding level. In contrast to what 
happens in vacuum, the screening of the conduction elec- 
trons reduces the bonding-antibonding energy splitting, 
and the antibonding level may be occupied by conduction 
electrons if the Fermi energy of the metal is high enough 
[Fo|. The occupation of antibonding states weakens the 
H 2 bond and increases the bond length. Our results agree 
qualitatively with the jellium model calculations which 
also found the H2 binding energy to be positive and the 
bond length increased, ranging from 0.81 to 0.86 A de- 
pending on the jellium density. In particular, for low 
jellium electron density (corresponding to the center of a 
vacancy site in Al), the binding energy was found to be 
+0.02 eV |l6|. Similarly, one can calculate the binding 
energy of multiple H atoms trapped at a single vacancy 
site, which turn out to be positive as well. Based on these 
results, we conclude that if the single vacancy concentra- 
tion Cy is greater than the H concentration Ch, each 
vacancy in equilibrium should contain no more than one 
H atom. 

On the other hand, if Ch is greater than Cy, the ques- 
tion arises as to where will the extra H atoms be situated, 
at interstitial or at vacancy sites? Experimental measure- 
ments for the ratio Ch/Cv in Al range from 0.25 to 4, 



depending on H charging conditions, with the most prob- 
able value close to 1 0. To answer the above question, 
we have calculated the trapping energy E tra p of multiple 
H atoms at a single vacancy site, which is defined as: 



trap{n) = -[E c (VAl 



-nR)-E c (V M )}-[E c (U)-E% 

(2) 

where E c (Vai + JiH) is the cohesive energy of a system 
with n H atoms each situated at a single vacancy site, 
E®(H) is the cohesive energy of bulk Al with a H atom 
at the tetrahedral interstitial site, and E® is the cohe- 
sive energy of the ideal bulk without H. A negative value 
for the trapping energy represents the energy gain when 
the H atoms are trapped at a single vacancy site relative 
to being dispersed at n different tetrahedral interstitial 
sites. The results for E tra p as a function of n are sum- 
marized in Fig. 2. Consistent with the binding energy 
calculations, it is energetically most favorable for each 
vacancy to trap a single H atom. At the same time, it 
is also energetically favorable for multiple H atoms to 
be trapped at a single vacancy site relative to being dis- 
persed at interstitial sites as individual atoms. In fact, 
up to twelve H atoms can be trapped at a single vacancy 
in Al, twice the highest number of H atoms (six) that can 
be trapped in Fe [5j. The atomic arrangement of the 12 
H atoms trapped at a single vacancy is indicated in Fig. 
1(b). There are two H2 units in each (100) direction sur- 
rounding the vacancy, and the bond length is 1 A for all 
six units. The inter-molecule distance in each direction 
is 3 A, and the NN distance between H and Al in each 
direction is 2 A(the lattice constant of Al is 3.99 A). The 
ordered arrangement of the H atoms is necessary to min- 
imize the electrostatic energy. The greater H-trapping 
capacity of Al compared to Fe, can be attributed to its 
larger lattice constant and more delocalized nature of 
electrons. It is observed that the volume change of the 
supercell owing to the H additions is negligible. 

The fact that the H2 units at a single vacancy site 
attract the conduction electrons from the edge of the 
vacancy, raises the interesting possibility that the cova- 
lent/metallic bonds between the first shell of NN Al ions 
around the vacancy site may be disrupted enough to per- 
mit a coalescence of multiple vacancies. To check this 
possibility, we carried out calculations for a number of 
relevant configurations. Specifically, we have examined: 
(i) two vacancies, each with one H atom, forming a NN 
divacancy with two H atoms trapped; (ii) n vacancies, 
each with two H atoms, forming a complex of NN multi- 
vacancies with 2n H atoms trapped, for n = 2 and 3. To 
summarize the results, we use the notation of chemical 
reactions: 

2 VaiH -> (Vai) 2 H 2 - 0.21 eV (i) 
nV Al U 2 -> (V A i)„H 2 „ + n 0.29 eV (ii) 

where the last number in each equation represents the 
reaction enthalpy AH. A positive value of AH means the 
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reaction is exothermic, that is, the process from left to 
right is energetically favorable. AH is defined as follows 
for reaction (i): 

AH = 2E C {V A1 R) - £ c [(y A i) 2 H 2 ] - E° c , (3) 

where ^[(Vai^Hs] is the cohesive energy of a system 
with two H atoms trapped at a divacancy, with analogous 
definitions for reaction (ii). Consistent with our earlier 
discussion, we find reaction (i) to be unfavorable (en- 
dothermic) because the effect of a single H atom on the 
covalent/metallic bonding of the NN Al atoms around the 
vacancy site is small and localized. On the other hand, 
reaction (ii) is favorable for n = 2 and 3, because the 
H2 units can attract more conduction electrons from the 
nearby Al atoms, weakening the bonding among the NN 
Al atoms, which in turn drives the formation of multi- 
vacancies. The large energy gain in forming the tri va- 
cancy (n = 3) is of particular interest. First, it is con- 
sistent with the experimental observation that the single 
vacancy defects occupied by H atoms can coalescence to 
form platelets on {111} planes of Al. Although our calcu- 
lations primarily concern the formation of the trivacancy, 
it is likely that even larger vacancy clusters can also be 
formed based on the same mechanism. In support of this 
claim, we mention that the increase in positive enthalpy 
associated with reaction (ii) is linear in the number of va- 
cancies for n = 2 and 3. Second, these vacancy clusters 
can serve as embryos of cracks and microvoids with lo- 
cal H concentrations much higher than the average bulk 
value. 

Next we discuss the implications of our results on hy- 
drogen embrittlement phenomena. It was generally be- 
lieved that H-induced embrittlement in metals takes the 
form of plastic rupture rather than brittle fracture, con- 
sistent with the notion of H-enhanced local plasticity 
(HELP). It was widely observed that the fracture surface 
is along the active slip planes where shear localization oc- 
curs. For fee metals, the slip planes are the {111} planes. 
In many cases, microvoids open up along these active slip 
planes in front of the crack tip; these microvoids can open 
and close in response to the local stress. Plastic rupture 
occurs when these microcracks are joined to the crack tip, 
upon reaching the critical stress. Our results clearly sug- 
gest that H-enriched microvoids may be created along the 
slip planes by the coalescence of vacancies with trapped 
H. These microvoids can be formed only in the pres- 
ence of H, which produces an additional source of mi- 
crocracks necessary for the H embrittlement. Moreover, 
the H-induced vacancy formation also facilitates dislo- 
cation climb, leaving behind vacancy rows in the highly 
deformed regions, which may contribute to the formation 
of microcracks as well. Our studies, taken together with 
the observed vacancy-enhanced dislocation glide 0,0], 
suggest that vacancies are also responsible for the HELP 
phenomena that are a prelude to H embrittlement 0. 
The fact that there is a strong binding between H and 



dislocation cores, and H can enhance dislocation motion 
along the slip planes > provide a means of rapid trans- 
port of H atoms to the crack front. On the other hand, 
the apparent lattice mobility of H atoms is also enhanced 
since multiple H atoms may be trapped at a single va- 
cancy. All these vacancy-based mechanisms contribute 
to the H embrittlement as they increase the rate of crack 
growth. Finally, the significant H trapping at vacancies 
provides a scenario by which drastic increase of local H 
concentration may occur without improbable accumula- 
tion of H at bulk interstitial sites 5j. This new feature 
resolves the long-standing problem of how a sufficiently 
high H concentration can be realized to successfully in- 
duce H embrittlement in materials such as Al, where the 
equilibrium H concentration in bulk is extremely low. 
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TABLE I: The vacancy formation energy, AHy', the binding 
energy for the divacancy, A//| v =2Ai/£ - AH2V, where the 
last term is the formation energy of the divacancy. The total 
energy with a H atom occupying the octahedral interstitial 
site is set to zero, relative to which the total energy of a H 
atom occupying the tetrahedral interstitial site, Et, and the 
total energy of a H atom trapped at a single vacancy, Ey 
are defined. The last two columns are LDA and GGA results 
from other theoretical calculations. All energies are given 
in eV. The experimental values marked by an asterisk have 
been called into question due to incorrect interpretations on 
the experimental part, see ref. [IJ for details. 

LDA GGA Exp. [IP] LDA GGA 

AHy 0.66 0.54 0.67±0.03 0.70 [9] 0.54 [9] 

AHly -0.06 -0.07 0.2*, 0.3* -0.07 [9] -0.08 [9] 

E T -0.07 -0.05 [13] -0.13 [13] 

E v -0.47 -0.46 [13] 



[20] Using the energy curve between equilibrium and 
metastable positions of H in bulk Al, we estimated the 
tunneling probability between equivalent equilibrium po- 
sitions to be 3 x 10 -8 , and thus quantum effects are not 
expected to be significant. Thermal effects (the config- 
urational and vibrational entropy) are well beyond the 
reach of first-principles calculations for this system. 



